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I.  DESCRIPTION  OF  THE  PROBLEM 

The  offensive  force  consists  of  A  identical,  perfect  attackers  which  we  think  of  as 
RVs.  The  defense  is  a  two-layer  defense  consisting  of  an  area  defense  and  terminal 
defense.  All  defensive  weapons  are  perfect,  so  they  destroy  an  incoming  attacker  if  aimed 
at  it.  The  targets  have  value  v(i)  for  i=l,...,T,  and  we  order  the  targets  so  that 

v(l)  >  v(2)  >  ...  >  v(T). 

The  terminal  defense  is  preallocated,  with  the  number  of  interceptors  at  target  i 
equal  to  v(i).  The  area  defense  can  operate  in  two  different  modes  -  (1)  without  impact 
point  prediction  and  (2)  with  impact  point  prediction.  In  the  first  case,  no  ipp,  the  area 
defense  selects  B  of  the  A  incoming  RVs  at  random  and  fires  interceptors  at  those  RVs  to 
kill  them.  (B  is  the  number  of  area  interceptors.)  The  RVs  surviving  continue  to  their 
respective  targets,  where  they  encounter  the  terminal  defense.  The  terminal  defense 
engages  the  RVs  and  if  any  RVs  penetrate,  the  target  is  destroyed.  In  the  second  case,  with 
ipp,  the  defense  can  see  the  structure  of  the  attack  and  then,  with  this  knowledge,  allocate 
the  area  interceptors  to  specific  RVs.  For  any  attack,  the  total  value  destroyed  will  be 
greater  in  case  (1)  than  in  case  (2),  assuming  that  the  defense  employs  the  best  strategy  in 
the  area  defense  layer.  The  problem  is  to  find  the  optimal  attacks  in  both  cases  (optimal 
meaning  greatest  expected  value  destroyed)  and  to  examine  quantitatively  the  difference 
between  the  two  cases. 

H.  A  SIMPLE  EXAMPLE  AS  A  GUIDE  TO  UNDERSTANDING 

In  order  to  illustrate  many  concepts  we  will  consider  the  following  fairly  simple 
situation: 

•  targets  of  value  10,  8,  5,  2, 1 

•  area  defense,  B=5 

•  number  of  attackers  variable. 


First  consider  the  case  with  no  ipp.  Suppose  the  attacker  has  10  weapons.  His 
problem  is  to  decide  how  to  distribute  the  10  weapons  so  as  to  maximize  the  expected 
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damage,  knowing  that  the  defense  does  not  have  impact  point  prediction.  Some  points  are 
clear: 

•  If  the  offense  elects  to  attack  target  number  i  then  it  is  useless 
to  aim  a  number  of  RVs  less  than  or  equal  to  v(i).  So  in  our 
case,  the  offense  should  not  attack  target  1  at  all,  if  he  attacks 
target  2  he  should  use  at  least  9  of  his  weapons,  etc. 

•  If  two  targets,  say  target  i  and  target  j  are  attacked,  then  if  i  <  j 
(so  that  v(i)  £  v(j)),  then  an  optimal  attack  must  have 

a(i)  -  v(i)  >  a(j)  -  v(j), 

where  a(k)  is  the  number  of  attackers  aimed  at  target  k. 

Professor  Soland  has  written  a  dynamic  programming  algorithm  to  determine  the 
optimal  targeting  for  the  attacker.  Appendix  A  contains  a  short  description  of  his 
algorithm. 

Using  his  algorithm  we  find  that  with  10  weapons  the  optimal  attack  is 

a(l)  =  0 
a(2)  =  0 
a(3)  =  0 
a(4)  =  7 
a(5)  =  3, 

and  the  expected  damage  is  2  1/3. 

I  have  programmed  the  algorithm  in  FORTRAN  on  IDA's  VAX  785  and  run  it  for 
different  numbers  of  attackers.  A  listing  of  the  program  and  the  outputs  from  its  various 
runs  can  be  found  in  Appendix  A.  The  cases  treated  there  are  of  area  defense  equal  5,  15, 
25,  35,  and  45. 

The  results  are  shown  in  Figure  1.  The  solid  curves  show  the  cases  B=5,  15,  25. 
In  each  case  the  shape  of  the  curve  is  similar,  with  a  staircase  at  low  attack  levels  and  a 
smoother  curve  sloping  to  the  right  at  higher  attack  levels. 

Consider  now  the  case  where  the  area  defense  has  impact  point  prediction.  When 
the  defense  has  ipp,  the  offense's  optimal  attack  changes.  There  are  two  basic  strategies 
that  I  have  defined  for  the  offense  (1)  shoot-to-kill  and  (2)  defense  dilution. 
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Figure  1.  Difference  in  the  number  of  attackers  needed  to  obtain  a  given 
damage  level  with  and  without  impact  point  prediction  in  the 
case  of  five  targets  valued  10,8,5,2,1. 


With  the  first  strategy,  the  offense  chooses  a  list  of  targets  to  attack  and  then  attacks 
a  target  of  value  v  with  v+B+1  RVs.  This  shoot-to-kill  attack  effectively  renders  the  area 
defense  impotent  to  defend  the  target  because  even  if  the  defense  were  to  use  all  of  his  B 
area  interceptors  to  defend  the  target,  v+1  RVs  will  still  penetrate  to  the  terminal  defense 
which  consists  of  v  interceptors  --  the  target  is  destroyed.  Thus,  at  each  attack  level,  the 
optimal  shoot-to-kill  strategy  consists  of  picking  the  largest  aggregate  target  value  to  attack 
subject  to  the  condition  of  having  enough  RVs  to  attack  a  target  of  value  v  with  v+B+1 
RVs. 

The  maximum  damage  possible  at  each  attack  level  for  B=5,  15,  25,  and  35  using 
the  shoot-to-kill  strategy  can  be  found  in  Appendix  B. 

It  turns  out  that  the  shoot-to-kill  strategy  is  optimal  for  the  offense  when  either  the 
attack  level  is  small  and/or  the  number  of  area  interceptors  is  small.  In  fact  for  B=5  and  15 
discussed  above,  the  shoot-to-kill  strategy  is  optimal  at  all  attack  levels. 

However,  when  B=25  (so  we  no  longer  have  B  "small")  and  A=62  (so  the  attack 
level  is  not  "small")  the  shoot-to-kill  strategy  is  not  optimal.  To  prove  it  is  not  optimal,  I 
shall  demonstrate  a  better  strategy.  The  optimal  shoot-to-kill  strategy  would  be  to  attack 
just  the  target  of  value  10  because  no  other  feasible  combination  has  a  higher  value 
destroyed.  However  consider  the  attack 

a(l)  =  10+13  =  23 

a(2)  =  8  +  13  =  21 

a(3)  =  5+  13  =  18 

a(4)=  0 

a(5)=  0. 

The  defense,  with  25  area  interceptors,  can  save  only  one  of  the  three  targets  attacked.  To 
minimize  his  losses,  he  will  save  the  value  10  target,  so  the  offense  destroys  the  targets  of 
value  8  and  5  for  a  total  damage  of  13. 

The  offense  has  diluted  the  defense  by  attacking  more  targets  than  with  the  shoot- 
to-kill  strategy.  The  offense  does  not  dictate  which  targets  will  be  destroyed  with  certainty 
but  leaves  that  choice  to  the  defense.  What  the  offense  can  guarantee  is  the  total  damage. 
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Further  properties  of  the  defense  dilution  strategy  will  be  examined  n  the  next 
section. 

By  examining  the  possible  attacks  for  the  cases  B=5, 15, 25  we  can  find  the 
optimal  attacks,  as  well  as  the  expected  damage.  These  results  are  plotted  on  Figure  1  as 
the  dashed  lines.  Note  that  for  a  given  number  of  area  interceptors,  the  curve  for  ipp  is 
always  to  the  right  and  below  the  curve  for  no  ipp. 

The  quantitative  comparisons  we  will  be  making  between  the  two  curves  will  take 
the  form: 

To  achieve  a  specified  damage  level,  what  percentage  more  RVs 

must  the  offense  expend  when  the  defense  has  impact  point  prediction,  as 

compared  to  the  case  when  the  defense  does  not  have  ipp? 

For  the  cases  B=5,  15,  25,  and  v=13,  20  the  results  are  shown  in  the  lower  right- 
hand  comer  of  Figure  1.  For  example,  with  15  area  interceptors  the  offense  needs  54.5% 
more  RVs  to  achieve  a  damage  level  of  20  when  the  defense  has  ipp  than  when  it  doesn't. 
Stated  another  way,  if  the  offense  can  deny  the  defense  impact  point  prediction  (for 
example  by  maneuvering  after  the  area  defense  layer)  then  the  offense,  to  achieve  a  damage 
level  of  20,  can  reduce  his  forces  by  approximately  1/3  (=  .54/1.54). 

HI.  THE  DEFENSE  DILUTION  STRATEGY 

When  the  defense  has  many  interceptors,  the  shoot-to-kill  strategy  is  nonoptimal 
because  the  cost  to  attack  a  new  target  is  high  relative  to  its  value.  What  is  called  for  is  a 
new  strategy  which  we  name  the  defense  dilution  strategy. 

Suppose  that  the  offense  attacks  all  of  the  targets  in  the  target  set,  meaning  that  a(i) 
>  v(i)  for  all  i.  The  area  defense  in  his  choice  of  which  targets  to  defend  (there  is  ipp)  must 
solve  a  knapsack  problem: 

^zj(a(i/-v(i)  <  B 


MAX 

i  £  {0,1} 


This  is  generally  a  difficult  problem  to  solve  exactly.  One  heuristic  solution  would  be  to 
order  the  "densities" 


v(i) 

a(i)— v(i) 


in  increasing  order,  and  then  let 


1  i^k 
0  i>k 


where 

x 

k  =  MAX  s.t.  Y  (a(i)-v(i))  £  B. 

*  i=l 

x  k 

The  value  saved  is  the^v(i).  However,  this  leaves  B  -^(a(i)-v(i))  area  interceptors 
unutilized,  so  that  we  suspect  this  is  not  optimal  for  the  defense. 

In  fact,  it  is  easy  to  make  up  examples  where  this  strategy  is  not  optimal. 
However,  this  strategy  results  in  a  value  destroyed  which  is  at  most  v(k+l)  more  than 
optimal,  so  if  the  target  values  are  small  relative  to  the  total  value,  then  the  percentage 
difference  between  the  two  is  small. 

We  can  consider  the  continuous  (or  partitionable  item)  knapsack  problem 


MAX  TzivO) 
zi£[0,l] 

I*i  (a(i)— v(i)>  <>  B 


which  has  the  solution 


zi  = 


1 

k 

B  -£(a(i)-  v(i)) 

i=l _ 

a(k+l)  -  v(k+l) 


i  <  k 


i  =  k+1 


\ 


The  answer  to  this  problem  is  at  most  v(k+l)  greater  than  the  original  knapsack 
problem  value.  So  again,  if  v(k+l)  is  small  relative  to  the  total  value,  the  percentage 
difference  is  small. 

The  offense's  problem  is  what  might  be  called  the  "knapsack  game": 

MIN  MAX 

5>G)=A  zi£{0,l} 

aO)>vO)  X  3  (a(i)-v(i))  £  B 


If  the  inside  problem  were  the  partitionable  item  knapsack  problem  (z,  £{0,1} 

instead  of  z,  e{0,l}),  which  we  assume  is  a  good  approximation,  then  the  solution  of  the 
outside  problem  is  to  take  aG)  such  that 


(Note  that  we  are  letting  a(i)  be  rational,  but  again  this  game  is  a  good  approximation  to  the 
game  where  a(i)  is  integral.) 


This  is  what  might  be  called  the , 


When  the  offense 


uses  this  attack,  the  value  of  the  continuous  game  is 


VALUE  _ 
SAVED  " 


Xv(i) 


B 

A-Iv(i) 


If  we  let  V  =  Xv(i),  then  we  have 


VALUE  _  VYl-  —  1 
DESTROYED  A-V' 


Note  however  that  at  low  attack  levels  the  offense  will  not  attack  all  of  the  targets  (the 
above  assumes  that  all  of  the  targets  are  attacked).  Let  us  assume  that  the  offense  will 
attack  a  fraction  f  of  the  total  value.  Then 


8(0  «  VALUE  _fV.fl— M- 
81  '  DESTROYED  V  A-fVJ 


The  offense  will  pick  f  so  that  this  is  maximized.  Taking  the  derivative, 

so  that  at  the  maximum 


fv  =  VA"  (VX  -  VB")  • 
if  v<s  vx  Oa -JS) 


then 


VAT  I  IF 
DESTROYED 


(Va-  Vb)2 


These  formulas  constitute  the  defense  dilution  strategy.  For  a  more  exact  analysis  of  this 
knapsack  game  and  the  defense  dilution  strategy,  see  Appendix  C.  The  algorithm  there 
was  never  implemented  -  instead  it  was  used  as  justification  of  the  two  formulas  above  for 
the  expected  value  destroyed. 
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IV.  LARGE  SCALE  PROBLEMS  WITH  ONE  AREA 

In  the  process  of  solving  the  ipp  vs.  no  ipp  on  a  nationwide  scale,  we  will 
encounter  rather  large  size  problems,  in  the  number  of  targets,  number  of  interceptors  and 
number  of  attackers. 


The  range  of 


number  of  targets  is  43  -  80 
total  target  value  is  430  -  800 
number  area  int  is  43  -  800 
number  attackers  is  50  -  2000 


As  an  example  of  the  qualitative  results  for  large  scale  problems  compared  to  the  small 
example  in  section  n,  I  present  the  following  situation 


targets 


5  targets  of  value  20 
10  targets  of  value  10 
12  targets  of  value  5 
Total  value  =  260 

Number  of  area  interceptors  =  50,150,250. 


The  results  are  shown  on  Figure  2.  As  we  see,  the  no  ipp  curves  have  become 
nearly  lines  and  the  ipp  curves  have  smoothed  quite  a  bit,  except  at  the  lower  end  where  the 
shoot-to-kill  strategy  with  its  characteristic  staircase  shape  is  present. 

The  solid  curves  were  calculated  at  25  step  intervals  using  Soland's  algorithm  and 
the  dashed  curves  were  calculated  as  a  combination  of  stk  and  dd  strategies. 

I  cannot  account  for  the  linear  shape  of  the  no  ipp  curves  at  present.  However,  they 
seem  to  be  approximated  by  lines  joining  the  following  two  points 

(B,0)  and  ((B+T+V)  (1+a),  V) 


where  T  =  number  of  targets  and  a  is  some  number  between  .  10  and  .20. 


VALUE  OF  IPP 

TARGETS  =  20,  20,  20,  20,  20 


Figure  2.  Difference  between  impact  point  prediction  and  no  impact  point 
prediction  for  a  single  citadel. 


B  ■  1 1 M  I  ft  I  yj  liTSjHIlTO 


Tilt] 


MjE 


a{i>  (1--^  )  Q  V(i)  (IF  a(i)  a  v(i)) 


See  Appendix  F  for  a  little  further  explanation. 


Computationally,  it  would  be  a  great  benefit  to  have  some  good  approximation 
because  computing  the  no  ipp  curves  is  the  largest  computational  burden.  The  calculation 
of  one  of  the  curves  for  a  large  problem  was  taking  >  5  hours  of  CPU  time. 


Appendix  F  contains  a  thorough  analysis  of  the  structure  of  optimal  attacks,  which 
leads  to  an  approximation  for  the  expected  damage  in  the  no  ipp  case. 


V.  LARGE  SCALE  PROBLEMS  WITH  MANY  AREAS  -  NATIONWIDE 
DEFENSE 


We  have  30  nonoverlapping  areas  (which  we  call  citadels)  with  targets  in  each  and 
independent  area  defenses  in  each. 


The  target  structure  of  the  citadels  is  shown  on  Figure  3.  In  a  given  row  of  the 
table,  we  read  across  the  number  of  the  citadel(s)  with  the  particular  target  structure,  the 
number  of  targets  of  value  200,100,20,10,5  and  the  total  target  value. 


The  total  target  distribution  is 


5  targets  of  value  200 


10  targets  of  value  100 


140  targets  of  value  20 


800  targets  of  value  10 


1000  targets  of  value  5 


*/  "  •  V  '«*  *  .*  V  “  __  .*  •  .  '  .  •  .  *  ,  '  ^  1  1.^”  •/ 


for  a  total  of 


1955  targets  with  total  value  17,800. 

I  ran  three  different  cases  of  the  problem: 

A.  Area  defense  =  10%  of  target  value 

B.  Area  defense  =  40%  of  target  value 

C.  Area  defense  =  100%  of  target  value. 

In  each  case  I  computed  the  no  ipp  and  ipp  curves  for  each  of  the  individual  citadels 
and  then  used  dynamic  programming  with  a  step  size  of  50  to  determine  the  optimal  attack 
against  the  nationwide  target  set. 


Citadel  Number 

m 

Number  of  Targets  of  Given  Value 

1QQ  2Q  10  5 

Value 

1-6 

0 

0 

2 

30 

18 

430 

7-12 

0 

0 

5 

20 

50 

550 

13-18 

0 

0 

8 

30 

36 

640 

19-20 

0 

0 

0 

50 

30 

650 

21-23 

0 

1 

7 

20 

20 

540 

24-25 

0 

2 

3 

40 

28 

800 

26 

1 

0 

2 

20 

50 

690 

27 

1 

0 

3 

15 

40 

610 

28 

1 

0 

10 

15 

30 

700 

29 

1 

1 

4 

15 

50 

780 

30 

1 

2 

4 

15 

30 

780 

Totals: 

5 

10 

140 

800 

1000 

1955  Targets 

_ 17,800  Value _ 

Figure  3.  DATA  SET  USED  IN  THE  MULTI-CITADEL  CASE 
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The  results  are  shown  on  the  following  four  pages  as  Figure  4  and  Tables  1  and  2. 
There  is  approximately  a  33%  difference  in  the  no  ipp  and  ipp  curves.  Appendix  D 
contains  the  curves  comparing  no  ipp  and  ipp  for  the  individual  citadels. 


The  results  are  shown  in  the  same  format  as  those  of  the  10%  case,  and  are 
displayed  as  Figure  5  and  Tables  3  and  4.  There  is  approximately  a  75%-85%  difference 
between  the  no  ipp  and  ipp  curves.  See  Appendix  E  for  further  results. 


For  this  case,  because  of  the  lengthy  computation  time,  I  used  the  approximation 
derived  in  Appendix  F  for  the  no  ipp  curves.  There  is  a  125%  difference  between  the  np 
ipp  and  ipp  curves.  See  Appendix  F  for  the  results  in  this  case. 


20,000 


e  target  value. 


Table  1. 

Area  Defense 

=  10%  of  Value 

Value 

Number  of  Attackers 
no  ipp  ipp 

Difference 

Percentage 

Difference 

4,450 

(1/4) 

5,900 

7,900 

2,000 

33.9% 

8,900 

(1/2) 

12,150 

16,200 

4,050 

33.3% 

13,350 

(3/4) 

18,450 

24,750 

6,300 

34.1% 

17,800 

(all) 

27,450 

The  percentage  difference  between  the  number  of  attackers  needed  with  and  without  ipp  to 
obtain  three  different  damage  levels  when  the  area  defense  equals  10  percent  of  the  target 


Table  2.  TARGET  SET  OF  30  CITADELS 
Area  Defense  =  10  percent  value 


Attackers 

no  ipp 

lEE 

1000 

773.67 

630.00 

2000 

1537.42 

1249.74 

3000 

2284.42 

1817.14 

4000 

3031.42 

2351.64 

5000 

3777.03 

2888.61 

6000 

4515.46 

3425.91 

7000 

5251.23 

3962.14 

8000 

5976.14 

4498.07 

9000 

6684.60 

5034.15 

10000 

7393.58 

5570.92 

11000 

8102.04 

6106.84 

12000 

8810.53 

6642.76 

13000 

9516.66 

7178.31 

14000 

10223.45 

7714.23 

15000 

10926.22 

8251.00 

16000 

11633.89 

8785.77 

17000 

12336.13 

9322.53 

18000 

13034.26 

9857.24 

19000 

13728.13 

10392.85 

20000 

14420.93 

10928.71 

21000 

15118.95 

11463.16 

22000 

15807.29 

11995.89 

23000 

16500.12 

12508.39 

24000 

17131.27 

12986.21 

25000 

17534.75 

13472.32 

26000 

17721.97 

13934.17 

27000 

17788.48 

14347.01 

28000 

17800.00 

14686.42 

29000 

17800.00 

14968.50 

30000 

17800.00 

15195.92 

31000 

17800.00 

15397.85 

32000 

17800.00 

15565.93 

33000 

17800.00 

15714.10 

34000 

17800.00 

15841.92 

35000 

17800.00 

15956.75 

Continued 
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TABLE  2.  CONCLUDED 


Attackers 

36000 

37000 

38000 

39000 

40000 

41000 

42000 

43000 

44000 

45000 

46000 

47000 

48000 

49000 

50000 


noipp 

17800.00 

17800.00 

17800.00 

17800.00 

17800.00 

17800.00 

17800.00 

17800.00 

17800.00 

17800.00 

17800.00 

17800.00 

17800.00 

17800.00 

17800.00 


ipp 

16057.79 

16148.92 
16230.27 
16305.14 

16372.81 

16433.81 
16490.35 
16542.21 
16590.32 

16634.82 

16676.11 

16714.71 

16750.64 

16784.30 

16815.93 


Expected  damage  levels  for  different  numbers  of  attackers  with  optimal  attacks  and 
defenses  in  the  cases  of  ipp  when  the  area  defense  equals  10  percent  of  the  target  value. 


Vs  s  s 


I 


t' 
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target  set  when  the  area  defense  is  equal  to  40  percent  of  the  target 


Table  3.  Area  Defense  =  40%  of  Value 


Value 

Number  of  Attackers 
no  ipp  ipp 

Difference 

Percentage 

Difference 

4,450 

(1/4) 

7,850 

13,900 

6,050 

77.1% 

8,900 

(1/2) 

16,000 

28,600 

12,600 

78.8% 

13,350 

(3/4) 

24,400 

45,000 

20,600 

84.4% 

17,800 

(all) 

35,150 

The  percentage  difference  between  the  number  of  attackers  needed  with 
and  without  ipp  to  obtain  three  different  damage  levels  when  the  area 
defense  equals  40  percent  of  the  target  value. 


TABLE  4.  TARGET  SET  OF  30  CITADELS 
Area  Defense  =  40  percent  of  value 


Attackers 

no  ipp 

ipp 

1000 

571.06 

400.00 

2000 

1141.96 

800.00 

3000 

1708.47 

1100.25 

4000 

2277.82 

1438.75 

5000 

2852.13 

1740.99 

6000 

3414.18 

2046.36 

7000 

3977.21 

2349.46 

8000 

4532.07 

2652.86 

9000 

5087.93 

2956.15 

10000 

5646.87 

3259.54 

11000 

6188.95 

3562.59 

12000 

6744.81 

3865.99 

13000 

7277.76 

4169.06 

14000 

7817.72 

4472.23 

15000 

8352.77 

4775.65 

16000 

8885.43 

5078.92 

17000 

9424.26 

5382.19 

18000 

9963.09 

5685.22 

19000 

10491.94 

5988.34 

20000 

11026.02 

6291.61 

21000 

11564.85 

6594.88 

22000 

12091.34 

6897.80 

23000 

12623.87 

7201.03 

24000 

13154.54 

7504.30 

25000 

13682.89 

7807.57 

26000 

14205.98 

8110.49 

27000 

14729.07 

8413.49 

28000 

15252.16 

8716.54 

29000 

15775.25 

9019.54 

30000 

16298.34 

9322.54 

31000 

16821.43 

9625.44 

32000 

17251.72 

9928.51 

33000 

17251.72 

1023.51 

34000 

17692.32 

10529.98 

35000 

17791.93 

10830.37 

Table  4.  CONCLUDED 


Attackers 

no  ipp 

m 

36000 

17800.00 

11127.66 

37000 

17800.00 

11418.10 

38000 

17800.00 

11705.17 

39000 

17800.00 

12002.46 

40000 

17800.00 

12251.41 

41000 

17800.00 

12519.97 

42000 

17800.00 

12759.36 

43000 

17800.00 

12973.59 

44000 

17800.00 

13169.69 

45000 

17800.00 

13348.16 

46000 

17800.00 

13511.36 

47000 

17800.00 

13662.50 

48000 

17800.00 

13801.10 

49000 

17800.00 

13931.17 

50000 

17800.00 

14051.49 

Expected  damage  level  for  different  numbers  of  attackers  with  optimal  attacks  and 
defenses  in  the  cases  of  ipp  and  no  ipp  when  the  area  defense  equals  40  percent  of  the 
target  value. 


VI.  COMPUTER  PROGRAM  DOCUMENTATION 


All  of  the  calculations  were  done  on  the  VAX  785  system  at  EDA  with  programs 
written  in  FORTRAN. 

The  no  ipp  programs  were  run  in  a  batch  mode.  Since  I  didn't  (and  still  don’t) 
know  how  to  read  data  from  a  File  during  the  execution  of  a  FORTRAN  program,  each 
different  batch  job  consists  of  the  same  basic  program  that  does  the  work,  with  different 
"front  ends"  contouring  the  data  used  in  that  particular  job.  The  .for  program  is  explained 
in  Appendix  A.  The  other  files  are  self  explanatory. 


In  case  other  people  would  like  to  access  my  files,  here  is  a  diagram  of  my  directory 
structures: 


In  the  subdirectory  [.DEFENSE],  the  following  files  may  be  found: 
DYNAM10.COM  -  compiles,  links  and  runs  DYNAMIO.FOR 
DYNAM10.FOR  -  uses  dynamic  program  to  integrate  all  30  citadels 
THE10ANSWER.DAT  -  output  from  DYNAM10.FOR  attack  step  size  is  50 
THESHORT10ANSWER.DAT  -  shortened  form  of  the  above  file  -  step  size  =  1000 


DYNAM40.COM 

DYNAM40.FOR 

THE40ANSWER.DAT 

THESHORT40ANSWER.DAT 


same  as  above  except  now 

area  defense  =  40%  of  value  instead 
of  10%  of  value 
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In  the  subdirectory  [.BATCH.PER10],  we  find  the  files 

RUN1QNUMB%%%.* 

where  %%%  is  001,  002,  003, ...,  01 1  and  *  is  .com,  .dat,  .for,  .log. 

The  numbers  %%%  refer  to  the  data  sets  from  the  list  of  citadel  structures  (there  are 
1 1  citadel  structures  and  1 1  different  file  types).  Also,  area  defense  =10  percent  of  value 
for  each  citadel. 

Similarly,  the  subdirectories  [.BATCH.PER40]  and  [.BATCH.PER100]  contain 
the  files 

RUN4QNUMB%%%.*  and  RUN100NIJM%%%.* 

0 

where  for  the  first  of  those,  area  defense  =  40  percent  of  value,  and  for  the  second,  100 
percent  of  value. 


The  following  is  a  short  description  of  a  dynamic  programming  algorithm  due  to 
Prof.  Soland  that  appeared  in  an  IDA  memo  dated  18  July  1985.  Because  we  are 
concerned  with  the  case  of  perfect  attackers  and  perfect  defenders  there  are  some 
simplifications. 


There  are  A  attackers,  B  area  defenders,  T  targets  with  the  value  of  target  i  equal  to 
v(i),  and  at  target  i  there  are  v(i)  terminal  interceptors. 


For  i=l,2,...,T  and  x=0,l,...,A,  let 


W-  (x)  =  J  expected  target  value  destroyed  at  target  1 
* ' if  x  RVs  are  assigned  to  it  by  the  attacker.  J 


Since  at  target  i  there  are  v(i)  perfect  terminal  defenders,  we  see  that 


I  the  probability  that  more  than  v(i)  RVs  arrive  at  the] 
W;  (x)=  v(i)  •  ■%  terminal  defense  of  target  i,  given  that  x  RVs  were 

1  V  7  I  • _ j  .  •.  I 


,  aimed  at  it 


An  elementary  "ball  and  urn"  argument  shows  that  the  distribution  of  the  RVs  through  the 
area  defense  at  each  target  is  hypergeometric,  so 


x  (A-B)  (  B  ) 

•w.00  =  v(i)  I  y 

y=v(i)+l  fA\ 


To  find  the  optimal  attack  we  first  calculate  Wi(x)  for  i=l,2,...,T  and  x=0,l,...,A, 
and  then  apply  a  dynamic  programming  algorithm.  In  particular,  let 


s(u/)  = 


MAX  £wi(Zi) 
u  i=l 

£zi  =  r 
i=l 


Zj  a  non-negative  integer 


for  u=l,2,...,T  and  R=0,1,...,A.  We  can  calculate  these  s(u,r)'s  recursively  by  letting 


s(U)  =  Wi(r) 


for  r=0,l,...,A,  and  then 


v;- 

V.V-V 


s(u,r)  =  MAX  J  Wu(z)+  s(u-l,r-z)  s. 
z=0,l,...,r  \  J 

for  u=2,3,...,T  and  r=0,l,...,A. 

The  answer  we  are  seeking  is  then  just  s(T,A). 

This  algorithm  is  implemented  in  the  FORTRAN  program  reproduced  as  Figure  A- 
1.  It  is  the  basic  program  used  to  calculate  all  of  the  curves  in  the  later  appendices  when  the 
area  defense  does  not  have  impact  point  prediction . 


*r\.  <r 


DSA2:L’lFlNN.UNCLAS,i.dArCH.r>EA10JAUN10NwN3  001.FOii;i 

c 

r  DECLARING  THE  VARIABLES 

REAL  wfiAR(0:100,  0:5000 
SEAL  SUH(Q:10G,0: 5000) 

REAL  V (0: 1 00) 

REAL  S (0:5000) 

REAL  I,J,L,X,V  ,INIT,B,N,T 
REAL  A,C,D,K,PKAS£A,M00N,Z 
SEAL  LBIN/ANSW 


OPEN  <UNIT=10,  FIL£=‘SUN10NUM3001  .DAT*,  ST*  TUS=  *  NS  V  * ) 


2 

3 


4 


1 


C 

c 

c 


41 

42 


I N I T=  5C 
3  =  43 
N=50CC 
STEP=5C 
T  =  5  0 

DO  2  1=1,2 

V ( I )=20 

CONTINUE 
DO  3  1=3,32 

V  <I)  =  1 0 

CONTINUE 
DO  4  1=33,50 

VCI)=5 

CONTINUE 

HAXV=V<1) 

V(0)*V(1>+9  .  m 

y  RITE  (10,*)*THERE  ARE* , B, *  AREA  INTERCEPTORS  WITH  PK  *’,PXAREA 
WRITE  (1 0, *) *  T  HE  TARGET  VALUES  ARE* 

WRITE  (10,*)  <V(I),  1=1, T> 

WRITE  (10,*)  *  * 

DO  1  1=1, T 

VALUE=VALUE*V(I) 

CONTINUE 

DO  1C0C  A=INI T ,N, STEP 
ANSW=0 

THE  CALCULATION  OF  W9AR(I,J)  FOR  PERFECT  AREA  INTERCEPTORS 
(THERE  ARE  SHORTCUTS  IN  THIS  CASE  AS  COMPARED  TO  THE 
GENERAL  CASE.) 

DO  *2  X=0,  V (0 ) 

Off  41  J=0, A 
SUM ( X,4 ) =0. 0 
CONTINUE 
CONTINUE 
DO  90  1=1, T 

X=V(I) 

r=v«i-i) 

IF (r-X.LE. 10** (”3)  )THEN 
SOTO  90 

ENDIF 

DO  60  J=X*1,X+S 

DO  55  Z=X*1,NIN< J,T) 

IF(A-B-.5.LT.I)60T0  55 
IF(8-.5.LT.J-Z)60T0  55 


A-3 


0HA2:  L*F  INfc.  UNCLAS  >.d  ATCH  .P  =  H1  U.1W  U*1  JNUMJOOI  .F  OR;  1 
0RA2:  CNF  INN. UNCLASS. 3 ATCH. PE  9103 RUN10NUM8001. FORI  1 
CRA2:CMFINN. UN  CL  ASS. B ATCH. P E 9 1 Q] RUN1 ONUM 3 001 . FO A; 1 

IF (A-. 5.LT. J ) 60T0  55 

SUM<X,  J)=SUM(X,  J)*EXP(L3IN(  A-B,Z)  ALB  IN  (  3,  J -2  J-L3IN  (A,J  )  ) 
5  CONTINUE 

SUM<X,J>=SUN<X,J>*SUM<Y,J) 
IF<SUM(X,J).GE.0.999)G0T0  SO 
C  CONTINUE 

00  70  J=X*1*B,A 

SUM(X,J)=1 

C  CONTINUE 

GOTO  9C 

C  00  35  S=J,A 

SUM<X,3)=1 

5  CONTINUE 

0  CONTINUE 

00  92  1=1, T 

00  91  J  =0, A 

uaAR(t,j)=vm*sup(v(i),j) 

1  CONTINUE 

2  CONTINUE 


DYNAMIC  PROGRAMING  ALGORITHM  TO  FIND  THE  OPTIMAL  ATTACK 
00  120  9=0, A*. 5 

S(R)=W8AR <1,R) 

20  CONTINUE 

00  150  U*2, T- 1 ♦. 5 

00  140  R  =  0,A«-.5 
00  130  Z=0,A-R*.5 

$< A-R)*MAX(SC A-R),W3AR(U,Z) *S< A-R-Z) ) 
30  CONTINUE 

40  CONTINUE 

50  CONTINUc 

00  160  Z=0,A*.5 

ANSW*MAXC ANSW,W9ARCT,Z)*S< A-Z>) 

60  CONTINUE 


JO  00 
I  COO 


WRI TNG  THE  OUTPUT — EXPECTED  DAMAGE  8T  AN  OPTIMAL  ATTACK 

WRITE  <  UN IT*1 0,  FNT*2000)  A,  ANSW, <S«RT< A)-SORT< 9*1 ) ) • *2 
F0RMAKF1  2.0,  2  FI  2.2) 

IFCANSW. GE. VALUE*. 98JG0T0  9999 
CONTINUE 


•999 


END 


THIS  FUNCTION  CALCULATES  THE  LOGARITHM  OF  BINOMIAL 
COEFFICIENTS  3T  USING  STIRLING’S  APPROXIMATION 
TO  FACTORIALS 


REAL  FUNCTION  L3IN(X,Y) 
IF(X.LT.Q.5)THEN 


-_T.  t*  TV  TV  T\i  w- TV'-T.  v.n^ .-y'. 


DRA2:l"FINK.  UNCLASS.dArCN.PcRIuJRUNlONUNdijOl.fOR;! 
’DRA2:C'IFINN.UNCLASS.3ATCH.PER10JRUNl0NUMeC01.F0R;i 
"  D  R  A  2 :  C’lF  INN.UNCLA3S.3ATCH.PES103RUN10NUM3C01  .  f  or  » i 
L3IN=G.G 

ELS E  IF<Y.LT.0.5)TH£N 
L3 IN  =  0. 0 

ELSE  IF<X-V.LT.0.5)THSN 
LSI N =0.0 

cLSc 

L3Ih=L0G<SaRT <2*3.141 5326*X))+X*L0G<X)-X*<1/<12*X)> 

L J  IN =L BIN- <L0G<S8RT<  2*  3.1415926*Y  )) >Y *L0G ( Y)- Y*< 1 /<12*Y>>) 
L3IN=L3IN- <L0G<S3RT< 2 *3.1 415926 *<X-Y>>)) 

L3  IN=L3IN-<<X-Y)*LQG<X-Y)-<X- Y) *< 1/<12*<X-Y) )) > 

END  I  F 
RETURN 


The  following  presentations,  Table  A-l,  are  from  the  FORTRAN  program  in  the  case 
where  the  target  values  are  10,8,5,2,1.  This  particular  case  is  the  one  discussed  in  Section 
II  in  the  body  of  the  paper. 

TABLE  A-l 


Expected  Value  Destroyed  by  an  optimal  attack  when  the  area  defense  is  5  inter¬ 
ceptors  without  ipp. 


Area  Defense=5 


Numher  of  Attackers 


Expected  Value  Destroyed 


16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 


10.00 

10.00 

10.51 

10.85 

11.50 
11.87 

12.27 
13.07 
14.07 

14.80 
15.02 

16.81 

17.59 

17.91 
18.41 
19.12 

19.59 
20.33 

21.51 

22.27 

22.60 
23.18 

23.85 
24.36 

24.79 
25.22 
25.46 
25.66 

25.80 

75.91 


TABLE  A-l 


Expected  Value  Destroyed  by  an  optimal  attack  when  the  area  defense  is  15  inter¬ 
ceptors  without  ipp. 

Area  Defense  =  15 


26 

10.00 

27 

10.00 

28 

10.21 

29 

10.47 

30 

10.79 

31 

11.30 

32 

11.65 

33 

12.03 

34 

13.08 

35 

13.91 

36 

14.46 

37 

15.14 

38 

16.38 

39 

17.11 

40 

17.57 

41 

17.84 

42 

18.44 

43 

18.98 

44 

19.61 

45 

20.62 

46 

21.35 

47 

21.93 

48 

22.31 

49 

22.81 

50 

23.38 

51 

23.87 

52 

24.17 

53 

24.60 

54 

24.99 

55 

25.23 

56 

25.42 

57 

25.58 

58 

25.70 

59 

25.80 

60 

25.86 

Continued 
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TABLE  A-l 


Expected  Value  Destroyed  by  an 

optimal  attack  when  the  area  defense  is  25  inter- 

ceptors  without  ipp. 

Area  Defense  =  25 

Number  of  Attackers 

ExDected  Value  Destroyed 

36 

10.00 

37 

10.00 

38 

10.11 

39 

10.33 

40 

10.59 

41 

11.09 

42 

11.47 

43 

11.71 

44 

12.69 

45 

13.55 

46 

14.13 

47 

14.63 

48 

15.77 

49 

16.65 

50 

17.21 

51 

17.56 

52 

17.91 

53 

18.49 

54 

18.93 

55 

19.83 

56 

20.64 

57 

21.32 

58 

21.80 

59 

22.16 

60 

22.68 

61 

23.21 

62 

23.62 

Continued 


TABLE  A-l 


Number  of  Attackers  Expected  Value  Destroyed 


63 

23.95 

64 

24.39 

65 

24.70 

66 

24.98 

67 

25.20 

68 

25.39 

69 

25.53 

70 

25.65 

71 

25.74 

72 

25.81 

73 

25.85 

74 

25.89 

75 

25.92 

Continued 


TABLE  A-l 


Expected  Value  Destroyed  by  an  optimal  attack  when  the  area  defense  is  35  inter¬ 
ceptors  without  ipp. 
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46 

10.00 

47 

10.00 

48 

10.07 

49 

10.27 

50 

10.52 

51 

10.98 

52 

11.37 

53 

11.63 

54 

12.49 

55 

13.34 

56 

13.95 

57 

14.37 

58 

15.50 

59 

16.35 

60 

16.98 

61 

17.39 

62 

17.65 

63 

18.19 

64 

18.67 

65 

19.39 

66 

20.23 

67 

20.89 

68 

21.44 

69. 

21.88 

70 

22.23 

71 

22.78 

TABLE  A-l  Concluded 


Expected  Value  Destroyed  by  an  optimal  attack  when  the  area  defense  is  45  inter¬ 
ceptors  without  ipp. 

Area  Defense  =  45 

Number  of  Attackers  Expected  Value  Destroyed 

56 

10.00 

57 

10.00 

58 

10.05 

59 

10.25 

60 

10.48 

61 

10.92 

62 

11.30 

63 

11.57 

64 

12.36 

65 

13.19 

66 

13.82 

67 

14.26 

68 

15.29 

69 

16.18 

70 

16.82 

71 

17.26 

72 

17.56 

73 

17.99 

74 

18.49 

75 

19.09 

76 

19.93 

77 

20.64 

78 

21.21 

79 

21.66 

80 

22.02 

TABLE  A-l  Concluded 


Number  of  Attackers 

Expected  Value  Destroyed 

81 

22.50 

82 

22.97 

83 

23.39 

84 

23.72 

85 

24.08 

86 

24.41 

87 

24.69 

88 

24.93 

89 

25.15 

90 

25.31 

91 

25.41 

92 

25.57 

93 

25.66 

94 

25.73 

95 

25.79 

96 

25.84 

97 

25.87 

98 

25.90 

99 

25.93 

100 

25.94 

i 
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APPENDIX  B 

THE  SHOOT -TO-KILL  STRATEGY 


I  have  considered  the  following  situations  as  in  Section  II: 

*  target  values  10,  8,  5,  2,  1 

*  terminal  defense  of  10,  8,  5,  2, 1 

*  area  defense  of  5  with  impact  point  prediction 

*  attack  level  variable 

*  perfect  interceptors  and  perfect  weapons. 

The  problem  is  to  determine  the  optimal  attack,  i.e.,  the  one  that  has  the  highest 
expected  damage.  R.M.  Soland  has  written  an  algorithm  to  solve  this  problem.  There  is 
however  a  simple  procedure  to  produce  a  non-trivial  lower  hound  for  the  expected  damage. 

The  procedure  relies  upon  the  "shoot-to-kill"  as  "guaranteed  damage"  strategy. 
With  five  area  interceptors  available  to  the  defense,  the  offense  can  guarantee  damaging  the 
first  target  with  10  +  5  +  1  =  16  attackers  directed  at  it  Similarly,  the  second  target  will  be 
damaged  by  14  attackers,  the  third  by  1 1,  the  fourth  by  8,  and  the  fifth  by  7.  The  "shoot- 
to-kill"  strategy  is  the  one  where  the  offense  picks  the  targets  it  wishes  to  and  can  destroy 
with  its  available  attackers  and  then  applies  the  number  of  attackers  to  those  targets  in  order 
to  guarantee  killing  them.  Thus  the  problem  reduces  to  simply  determining  the  wisest 
choice  of  targets  to  attack. 

For  example,  with  20  attackers  the  offense  can  damage  the  following  combinations 
of  targets:  any  of  the  individual  targets;  targets  4  and  5;  targets  3  and  5;  targets  3  and  4. 
For  these  possibilities,  the  expected  damages  are  1,2,  5,  8,  10,  3,  6,  7.  So,  we  see  that 
the  best  attack  is  to  overwhelm  target  1,  for  an  expected  damage  of  10. 

Carrying  out  this  procedure  for  a  range  of  attackers  we  find  the  following: 
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A  TWO  LAYER  ONE  AREA  DEFENSE  WITH  IMPACT  POINT 

PREDICTION 


PERFECT  ATTACKERS  AND  PERFECT  DEFENDERS 


A.  Algorithm  for  the  Knapsack  Game  When  There  Are  Many  Objects  of 

Equal  Value 

We  shall  first  consider  the  case  where  all  of  the  targets  are  of  equal  value. 

There  are  t  targets  each  of  value  V.  There  is  a  terminal  defense  of  V  perfect 
interceptors  at  each  target  and  an  area  defense  of  B  perfect  interceptors.  The  area  defense 
has  impact  point  prediction,  which  means  that  the  area  defense  can  be  allocated  after  the 
defense  observes  the  attack. 

The  problem  we  will  solve  is  to  determine  the  minimum  number  of  RVs  the  attacker 
needs  in  order  to  guarantee  destroying  r  targets,  where  r  varies  between  one  and  L 

Mathematically,  the  problem  is  to  determine  the  least  A  such  that 


MAX 

y  a(i)  iS  A 
aft)  £N 


MIN  £f(a(i),d(i))  *  r- 

>i)  *  B 
i)e  N 


where 


f(a(i),d(i)) 


■dW)={o! 


1  if  V+d(i)  £  a(i) 
OifV+d(i)  <  a(i) 


N  =  {0,1,2,...}. 

It  is  easy  to  see  that  for  any  value  of  A,  among  the  set  of  optimal  attacks  will  be  one  having 
for  each  i 

either  a(i)  >  V  or  a(i)  =  0. 

Hence,  we  may  restrict  our  attention  to  such  attacks.  Suppose  then  that  the  attacker  attacks 
n  targets,  where  r  £  n  <,  l  If  the  attacker  is  to  guarantee  destroying  r  targets,  then  he  cannot 
allow  the  defense  to  save  (n-r)-f  1  targets.  Thus, 


•  »  *>  *  *  "  fc  '»“«*•  *  «  *  *  *  •  ’/  */  «*  **■  t",  •»’*.  •*,  m  .  •' 

v  %  V  %  •/  ’  JV  A  .*»  Aw  -V.%  A  i  -*•  A  m>  fA  *.  i:i 


(*) 


£(a(ia)-v)>B+l 

o 


for  any  n-r+1  distinct  subscripts  {icy}* 

If  we  add  the  (n-^+l^  different  inequalities  above,  we  find  that 

(**)  gXa(i)-v)  >  (B+1)- 

Let  us  renumber  the  targets  so  that 

a(i)  £  a(j) 

fori  <  j  . 

From  inequality  (*)  we  have 

(a(l)-v)  +  (a(2)-v)  +  ...  +  (a(n-r+l)-v)  >  B+l. 
By  the  pigeonhole  principle 

(a(n-r+l)-v)  a  . 

Thus, 

a(i)  -  v  £  for  i  >  n-r+1. 

In  order  to  minimize  the  required  number  of  attackers,  we  should  have 

n-r+1 

£  (a(i)-v)  =  B+l 
i=l 

with 


and 


a(i)  -  v  £  for  i=l,  ...,  n-r+1 

a(i)  =  f"iw+r|  for  1  =  n'r+2’  •••’ n- 

We  find  the  total  number  of  attackers  to  be 
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Attackers  =  nV  +  (B+l)  +  (r-1) 


To  get  the  final  answer,  we  minimize  this  function  over  n  e  (r,r+l, t}. 


It  would  be  convenient  to  have  an  analytic  approximation  to  the  minimum  number 
of  attackers  so  as  to  perform  analyses,  draw  graphs,  etc.  Therefore,  such  an 
approximation  will  be  derived. 

First  we  consider  the  continuous  game  -  where  a(i)  and  d(i)  are  non-negative  real 
numbers.  We  can  use  an  argument  similar  to  the  one  above  to  find  that  the  optimal  attack  is 
equally  distributed: 

a(i)  =  +  V  for  i  =  l,2,...,n 

and  so 

Attackers  =  n 

-  °V+  (B+l)  +  (r-1)  (2&). 

We  now  further  relax  the  constraint  of  n  being  integral.  Thus,  the  offense  is  free  to  choose 
the  exact  value  of  targets  to  be  attacked.  For  V  small  and  t  large  this  is  probably  a  valid 
assumption. 

In  any  event,  we  shall  use  the  following  combination  of  notations: 

Aint  refers  to  the  problem  where  a(i),  d(i)  e  N 
ACOnt  refers  to  the  continuous  problem 
A  refers  to  the  problem  where  n  e  N. 

Thus,  A[nt  would  be  the  answer  to  the  original  problem,  while  Acont  is  the  most  tractable 
analytically. 

To  find  Acont,  we  rewrite  the  number  of  attackers  as 

Attackers  =  (r-l)V  +  (B+l)  +  (n-r+l)V  +  (r-1)  =  f(n). 

The  Fust  two  terms  are  constant,  so  we  don’t  need  to  consider  them.  To  minimize 
the  latter  terms,  we  use  the  Arithmetic  Mean-Geometric  Mean  inequality  which  says  that  the 
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sum  of  two  positive  numbers  whose  product  is  a  constant  is  minimized  when  the  numbers 
are  equal,  which  implies 


unless  of  course 


(n*-r+l)V  .  (r-Djgil- 


n*-r+ 1  =  V(r-l)  (B+l)/V 


Acont  =(r"1)v  +  (B+l)  +  2V(r-l)  (B+1)V 
=  (V(r-l)V  +  V (B+l)  )2f 


n*-r+l  =  V(r-l)  (B+l)  /V  >  t-r+1 


in  which  case  the  minimum  occurs  at  n=t,  so 

Acont=  tV+ 

In  order  to  find  Acont,  we  know  that  n*  is  at  most  1/2  away  from  the  nearest 
integer.  And  we  know  that  f '(u)  is  decreasing,  so 

Acont  (f(n*+ l/2),f(n*-l/2))  =  f(n*-l/2). 

Substituting  this  in,  we  see  that 

Acont  ^  Acont  +  (f(n*-l/2)  -  A cont) 


If  n*  >  t,  then 


=  A  + 
cont 


1/4  V _ 

(r-1)  (B+l)  _  1/2 


=  A  +  - - - 

cont  4(n*-r+l/2) 


tf. ml — d 


Finally,  to  get  a  bound  on  Ajnt,  we  see  that  for  all  n, 


<t-»  (Srir)  *  [§&1  *  <r-»  (£m-)+  ('•!)• 


Therefore,  we  have 


Acont  -  Acont  -  Acont  +  ^r_1\ 


Combining  all  three  inequalities,  we  get 


Acont  S  Acom  £  Ata[  £  Acom  +  (r-l)<  Acom+  (r-1)  +  4(n*-r-l) 


and  if  n*  <  t, 


A  .  =  A  <A.  £  A  +  (r-1)  =  A.  +  (r-1). 

cont  cont-  mt  cont  int 


It  is  possible  to  derive  a  lower  bound  for  A .  by  considering  A . 


The  number  of  attackers  for  the  integral  solution  is 


Attackers  =  (r-l)V  +  (B+l)  +  (n-r+l)V  +  (r-1) 


which  we  wish  to  minimize  for  n  £  [r,t]. 


Letting  p  =  n-r+1,  so  that  p  e  [l,t-r+l],  the  number  of  attackers  can  be  rewritten 


Attackers  =  (r-l)V  +  (B+l)  +  pV  +  (r-l)j-E±i^ 


The  graph  of  this  function  is  a  series  of  line  segments  of  slope  V,  with  end  points  at 


P=  B+L  andP  =  ®±! 

1  l- 1 


■*  mJ*.  A.' 


The  minima  will  obviously  occur  at  c.+i  where  i  is  the  largest  integer  such  that 


c 


I  ^ 


v  (W  *  *  fr'1) 


i.e.. 


i*  = 


unless  ¥  >  t-r+1,  in  which 


case 


;* 


In  either  case, 


Aint=(i*+V)(r-1  +^±L) , 
which  is  another  lower  bound  for  A^t- 

B.  Algorithm  for  the  Knapsack  Game  When  There  Are  Many  Objects  of 
Unequal  Value 

Using  the  insight  gained  from  section  A,  we  can  now  outline  an  algorithm  for  the 
case  of  unequal  target  values. 

Suppose  the  targets  are  valued  v(i)  with 

v(l)  £  v(2)  >  ...>  v(T). 

One  way  to  solve  the  problem  is  to  postulate  an  attack  level  A  and  then  enumerate  all  of  the 
attack  vectors 

a  =  [a(l),  a(2), ... ,  a(T)J 

with 


0  <  a(i)  <  A 


and 


Xa(i)  =  A. 

For  each  one  of  these  attack  vectors,  the  defense  has  a  knapsack  problem  to  solve,  i.e., 

MAX  X  ZjVf 

such  that  Xzifi(a(i))  ^  B 

Zi£{0,l} 

where 

fi(x)  =  max(x-v(i),0). 

The  number  of  attack  vectors  to  be  examined  (and  hence  knapsack  problems  to  be 
solved)  is  on  the  order  of  A^'l,  which  is  beyond  feasibility  when  considering  attacks  on 
large  target  sets. 

In  the  case  of  perfect  attackers  and  interceptors,  it  is  simple  to  prove  that  optimal 
attacks  are  monotone  in  their  excesses  on  targets  attacked.  In  other  words,  if  i<j  (so  v(i)  > . 
v(j))  and  a(i)  =  0  a(j)  =  0,  then 

a(i)  -  v(i)  S  a(j)  -v(j). 

Using  this  we  can  diminish  the  number  of  attack  vectors  considered  by  roughly  a  factor  of 
T!,  so  the  number  of  knapsack  problems  to  be  solved  is  still  on  the  order  of  a(T-1). 

Sabbagh’s  implicit  enumeration  algorithm  does  not  fare  any  better  because  the 
"jumps"  in  the  implicit  enumeration  scheme  are  not  sufficiently  large  to  significantly 
decrease  the  number  of  attack  vectors  considered. 

The  algorithm  outlined  below  attempts  to  circumvent  the  computational  size  of  the 
above  algorithms  in  the  case  where  there  are  many  targets  of  the  same  value  by  means 
similar  to  the  ones  employed  in  part  A.  The  algorithm  must  solve  at  most 

G 

II  (n(i)  +  1) 
i=l 

linear  integer  programming  problems,  where  n(i)  is  the  number  of  targets  of  value  v(i)  and 
G  is  the  number  of  distinct  target  values. 
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I l- 


lc 


In  order  to  make  the  discussion  as  concrete  as  possible,  we  will  consider  the 
following  problem: 

6  targets  of  value  4 
4  targets  of  value  7 
Required  damage  level  =  25 
Area  defense  =  B 

We  wish  to  determine  the  minimum  number  of  attackers  needed  to  obtain  a  value  destroyed 

of  25.  | 

Figure  C-l  shows  the  set  up  for  the  algorithm.  Along  the  x-axis  is  the  number  of 
targets  of  value  7  and  along  the  y-axis  is  the  number  of  targets  of  value  4.  The  boundary 
line  for  u=25  satisfies  the  equation 


Points  to  its  left  represents  target  combinations  that  sum  to  less  than  25  value  and  points  to 
its  right  more  than  25.  The  19  integer  points  on  the  line  or  to  its  right  represent  feasible 
attack  points.  We  shall  examine  the  one  circled,  which  represents  the  attacker  attacking 
three  targets  of  value  7  and  4  targets  of  value  4.  Let  Xj(i=  1,2,3)  be  the  excesses  of  the 
attacks  against  the  targets  of  value  7  and  yj(j=l,...,4)  be  the  excesses  of  the  attacks  against 
the  targets  of  value  4. 

The  shaded  region  is  the  "defense  wins"  region,  which  means  that  if  the  defender 
can  save  a  combination  of  targets  to  get  within  the  region,  then  the  postulated  attack  does 
not  produce  a  destruction  of  25.  For  example,  if  the  defense  can  save  3  targets  of  value  4 
under  attack  and  2  targets  of  value  7  under  attack,  then  the  offense  destroyed  only  1  target 
of  value  4  and  1  of  value  7,  hence  a  total  of  1 1,  which  is  less  than  25. 

Hence,  points  in  the  region  represent  constraints  on  the  xj’s  and  yj's. 

The  minimal  points  represent  the  smallest  set  of  irredundant  constraints.  They  are 
defined  as  the  minimal  set  of  points  generating  the  whole  "defense  wins"  region,  where  the 
point  (a, (3)  "generates"  the  rectangular  set  of  points  (x,y)  satisfying  (x,y)  e  [0,a]  x  [0,(3]. 

The  uppermost  minimal  point  represents  the  constraint  that  the  defense  cannot  be 
allowed  to  save  any  two  of  the  targets  of  value  7.  Hence 
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C 


NUMBER  OF 
4  VALUE 

TARGETS  ATTACKED 


DEFENSE  WINS” 
EG 


Xj  +  x2£  (B+l) 

X!  +  x2£(B+1) 
x2  +  x3£(B+1). 

The  middle  minimal  point  represents  the  constraint  that  the  defense  cannot  be 
allowed  to  save  1  target  of  value  7  and  2  targets  of  value  4.  Hence 

xa  +  yp1  +  yp2  *  (B+1) 

for  a  e{  1,2,3}  and  Pi,P2  distinct  elements  of  { 1,2, 3 ,4]  . 

The  final  set  of  constraints  is  that  the  defense  cannot  be  allowed  to  save  all  4  targets 
of  value  4.  Thus, 

yi  +  y2  +  y3  +  y4  ^(B+i). 

What  we  need  to  do  now  is  minimize  the  total  attack, 

3  x  7  +  xj  +  x2  +  x3  +  4.x  4  +  yj  +  y2  +  y3  +  y4 

subject  to  these  constraints. 

By  summing  the  first  constraints,  dividing  by  2  and  adding  the  last  constraint  we 

find 

xi  +x2  +  x3+yi+y2  +  y3  +  y4>-^-(B+l) 

with  equality  if 

xj  =  1/2(B+1),  yj  =  1/4(B+1). 


All  of  the  middle  constraints  are  satisfied  with  these  values  of  x,  and  yj  so  we  have 
found  a  (potentially  non  integral)  optimal  attack.  If  B=3  (mod  4)  then  the  attack  is  integral. 
If  not,  we  can  get  good  bounds  on  the  number  of  attackers  required  by  decreasing  B  by  i 
and  increasing  B  by  4-i  where  B=i  (mod  4),  ie{  1,2,3}. 


Then  we  vary  the  feasible  attack  point  (circled  point)  to  one  of  the  18  other  feasible 
points.  At  each  point  we  solve  the  problem  and  take  the  minimum  of  all  the  answers  to  get 
the  final  result. 


This  algorithm  generalizes  immediately  to  more  than  2  distinct  target  values. 


C.  Remarks  and  Simplications  of  the  Algorithm 

1.  Suppose  the  feasible  attack  point  under  consideration  is  (a,p)  and  we  are  interested 
in  the  constraint  due  to  the  point  (a-a,  P~b).  Then  we  know  that 


a  (B+D 


where 


{i*}  is  any  collection  of  a  distinct  indices  in  { 1,2 . a}  and  {ja}  is  an  distinct  collection 

of  b  indices  in  {1,2, ...,  P}.  If  we  add  up  all  of  the  possible  combinations,  we  find  that 


(o-{0- -I-}*  (o  {(')•?  i’<Kxo 


(B+l) 


so  that 

2  (B+D. 

If  we  now  let  a,b,a,P  assume  continuous  values,  we  can  get  an  analytical  lower  bound. 
Suppose  we  take  a  and  b  so  that 

*-JL 

a  b 


and 


e  =  small  constant. 


(a-a)  x  7  +  (P-b)  x  4  =  25  -  e 


In  (he  limit,  a  -  a-  a„d  »  -  P"  7^4  ' 


7+4  p/a 
Then  the  constraint  becomes 


<Zxi+Syj)2  <b+u 


which  implies  ^xj  +  ^yj  - 


B+l 


1-  _2i~' 


la +  4p 


Thus 


Att  > 


B+l 


1  -  -25-,, 

7  a +4(3 

If  we  minimize  this  with  respect  to  a,P,  we  find  it's  a  minimum  when 


+  la  +  4  p. 


7a  +4p  =  25  V  25  (B+l) 


which  implies  that 


Att  > 


_E±L 


1  - 


_25_ 


■+  25  +  V  25  (B+l) 


25  +  V25  (B+l) 


=  (V2T  +  V B+T" )  , 

unless  25  +  V  25(B+1)  -  TOTAL  VALUE  =  4x7  +  6x4  =  52, 
in  which  case  we  would  take 

7a  +  4  p  =  TOTAL  VAL, 

so  that 


Att  > 


B+l 


1  + 


25 


+  TOT  VAL 


TOT  VAL 


=  TOT  VAL  J  1  + 


B+l 


{'  ’  TOTVAL  -  25 

=  52{,  + 


}• 


2.  The  algorithm  can  be  simplified  somewhat  if  the  exact  integral  answer  is  not  desired 
but  just  a  tight  bound  on  it. 


Suppose  there  are  N  distinct  target  values. 

We  shall  keep  the  same  terminology  as  in  the  case  of  two  target  groups,  except  we 

work  in  n-dismensional  space  where  the  boundary  line  becomes  a  boundary  hyperplane. 
Suppose  we  are  investigating  the  feasible  attack  point  (a1,a2,a3,...,an)  and  are 

concerned  with  the  constraints  due  to  the  vector  (bi,b2,...,bn).  As  above,  if  we  add  all  of 
the  possible  inequalities,  we  find 

Xj  C1  ~  Xx  i,  k  £  B+l. 

Let  C  be  the  convex  hull  of  the  lattice  points  in  the  "defense  wins"  region.  Let 

where 


Write  pas 


X  =  sup  {X:  X^xeC}  . 


n 

t  wj  •  (bj) 
j=l  J  J 

where  T?j  =  a  lattice  point  in  the  defense  wins  region 

OSwjS  l,£wj  -  !• 


i.e.,  find  p's  barycentric  coordinates. 

Since  each  of  the  bj's  is  a  constraint  vector,  we  have  for  all  j 


>  B+l. 


If  we  multiply  the  jth  equation  by  wj  and  add,  we  find 


1  '  oq  ^  Xxi,k  ^  B+1 


which  can  be  rewritten 


Thus  2<1-X)21*U*M 

and  ZSxi,k  ^  . 


Zl"bl’l  z2-bl,2  ...  zn'bl,n 
a2,l*bl,l  b2,2‘bl,2---b2,n-bl,n 


-D+A[Z^+  ..+AnZn=det 


bn,l'al,l  bn,2’al,2  bn,n’bl,n 
be  the  equation  of  the  hyperplane  through  the  points  b^. . . ,  bn. 


If  we  let 

£xi>k  =  -  _  B±i _ 

K  ’  (A1ai+..  +Ai/\xj+...  +An0Cn).  D 

Ai 

then  all  of  the  constraints  will  be  satisfied,  and  we  will  have  equality  in  the  constraint  for 
££xi>k  ,  so  we  have  obtained  the  minimum. 

However, 


C-15 


xi,k  =*4-  (B+l) 

where  t^sjeZ. 

This  is  the  rational  solution. 

The  integral  solution  must  have  Xj^eZ.  We  can  obtain  an  upper  bound  on  the 
number  of  attackers  by  letting 


C-16 


Then  the  total  number  of  attackers  assigned  to  the  ith  target  group  is 


Thus  we  have  the  bound 

“iVi+di  1l(B+1)  £  AttjS  |(aj- bj  jy  (B+l)l 

+  bi  pi-  <B+1>1 

<  cqvi  +  (ot i-bi)  ii-  (B+l)  +  1 
+  bi  iL  (B+i)  +  bj 
=  aivj+  (B+l)  +  (bi  +1) 


C-1 7 
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The  twelve  charts  contained  in  this  appendix  demonstrate  the  difference  between  ipp 
and  no  ipp  for  each  of  the  thirty  citadels  when  the  area  defense  equals  ten  percent  of  the 
target  value.  On  each  chart,  the  solid  line  represents  no  ipp  and  the  dashed  line  represents 
ipp.  The  citadels  covered  by  a  specific  chart  are  noted  at  the  top  of  the  chart. 


EXPECTED  VALUE  DESTROYED  EXPECTED  VALUE  DESTROYED 


EXPECTED 


The  eleven  charts  contained  in  this  appendix  demonstrate  the  difference  between  ipp 
and  no  ipp  for  each  of  the  thirty  citadels  when  the  area  defense  equals  forty  percent  of  the 
target  value.  On  each  chart,  the  solid  line  represents  no  ipp  and  the  dashed  line  represents 
ipp.  The  citadels  covered  by  a  specific  chart  are  noted  at  the  top  of  the  chart 


CITADELS  13-18 


CITADELS  21-23 


0  200  400  600  800  1000  1200  1400 


NUMBER  OF  ATTACKERS 


EXPECTED  VALUE  DESTROYED  EXPECTED  VALUE  DESTROYED 
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The  large  scale  problems  considered  in  this  paper  in  the  case  of  no  ipp  and  area 
defense  =  100%  of  value  presented  large  computational  burdens.  To  simplify  the 
computations,  I  derived  some  approximations  which  may  be  applied  with  benefit 


In  order  to  investigate  the  optimal  attacks,  I  supplemented  the  FORTRAN  program 
implementing  Soland's  algorithm  (discussed  in  Appendix  A)  by  the  introduction  of  the 
variables 


attop  (i,r)  =  number  of  attackers  used  against  target  i  in  an  optimal  attack  when 
the  offense  is  attacking  the  first  i  targets  with  a  total  of  r  RVs, 


numb  (i)  =  number  of  attackers  used  against  target  i  in  an  optimal  attack  against 
the  whole  target 


The  values  of  attop  (i,r)  are  obtained  readily  as  a  byproduct  of  the  dynamic  programming 
algorithm,  while  the  values  of  numb  (i)  can  be  easily  obtained  from  the  formulas 


numb  (T)  =  attop  (T,  a) 


for  i  =  T-l,  T-2,...,l. 


numb  (i)  =  attop  (i,  a  -  £numb  (j) ) 

j=itl 


After  examining  many  cases  of  the  optimal  attacks,  I  decided  upon  two  guidelines 
to  aid  in  the  approximations: 


1 .  If  the  attacker  attacks  two  targets  of  the  same  value  in  his  optimal  attack,  then 
he  attacks  them  with  equal  numbers  of  weapons. 


2.  When  possible,  the  attacker  attacks  more  valuable  targets  rather  than  less 
valuable  targets. 


Needless  to  say,  these  two  guidelines  held  true  in  the  many  cases  I  examined. 


At  this  point,  I  assumed  that  the  targets  were  equal  valued  and  numerous  enough 
that  the  attacker  is  not  limited  in  his  choice  of  targets.  The  arrival  probability  through  the 
area  defense  is  hypergeometric: 


nRVs  arrive,  /a_bN  /  b  N 

Prob  *  given  that  x  were  m  =  — n  Va\~  ' 
fired  at  the  target 


where  a  =  number  of  attackers  and  b  =  number  of  area  defenders.  I  approximate  this 
hypergeometric  distribution  by  a  binomial  distribution: 


n  RVs  arrive, 
Prob  -  given  that  x  were 
fired  at  the  target 


.apjBI2*.($)p*-nQn, 


where  p  =  b/a  and  q  =  1  -  b/a. 


Now  I  formulate  the  optimal  attack  problem  as 


MAX 


-t-X  (5)  (l-b/a)n  (b/a)*-" 


X  £{v+1,v+2,...}  n  =  r+1 

The  second  term  in  the  product  is  the  Incomplete  Beta  Function  for  which  I  found 
the  following  approximation: 

i \  „  a  +  b-l-n 

Ix(a,b)=£  (a  +  n_1)xn  (1-x) 
n  =  a 


z.  P(y)  +  e, 


where 


|£|  <  5  x  10"3  if  a  +  b  >  6  and  (a  +  b  -  1)  (1  -  x)  <,  .8 


3[wl(l-$b-w2  O-^a)] 


y  = 


(  wl2  +  ^2_)1/2 

~  a 

w1=(bx)1/3 
W9  =  (a  (1  —  x)) 

z  y 


1/3 


- 1/2 12 


dt 


Furthermore,  P(y)  has  the  following  rational  approximation: 

P(y)  =  1  -  1/2  (1  +  ajy  +  +  ^y3  +  a4y4)~4  +  e 

where 


|e|  <  2.5  x  1(H 


and 


aj  =  .196854 
a2  =  . 115194 
a 3  =  .000344 
34  =  .019527 


If  these  two  approximations  are  used  on  Iq  (v  +  1,  x  -  v),  then  we  can  numerically 
maximize  the  product 

-f  Iq  (v  +  1,  x  -  v) 

for  x  greater  than  v,  assuming  the  values  of  v  and  b/a  to  be  given. 

The  results  are  plotted  in  Figure  F-l  which  shows  p  vs.  b/a  for  different  values  of 
v.  The  quantity  P  is  defined  by 

6  =  100*  (  optimal RV s/target 

V  (v/q) 


F-3 
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Figure  F-2  shows  the  probability  of  kill  in  the  optimal  attack.  Together  these  two 
charts  are  the  basis  of  finding  the  approximate  optimal  attack.  As  an  example,  consider 
citadel  29  and  1400  attackers.  Using  Figure  F-l,  we  find  that 


Value 

Optimal  Number  of  Attackers 

4 

200 

503 

.9786 

100 

260 

.9685 

20 

58 

.9213 

5 

17 

.8286. 

With  1400  weapons,  the  attacker  will  use  503  against  the  200  value  target,  260 
against  the  100  value,  58  against  each  of  the  4  targets  of  value  20  and  31  against  13  targets 
of  value  10.  Two  RVs  will  be  left  over.  Summing  up  the  expected  damages  at  each  target, 
we  find  that  the  total  expected  damage  is  481.10. 

This  method  was  followed  to  generate  the  no  ipp  curves  for  each  of  the  individual 
citadels.  The  ipp  curves  were  generated  as  before  with  a  combination  of  shoot-tO-kill  and 
defense  dilution  strategies.  The  individual  citadels  were  then  integrated  into  a  nationwide 
system  by  using  dynamic  programming. 

Table  F-l  summarizes  the  differences  of  RVs  required  with  no  impact  point 
prediction  and  with  impact  point  prediction  in  selected  cases.  Table  F-2  gives  numerical 
details,  and  Figure  F-3  displays  the  aggregate  information.  The  remaining  figures  display 
resuits  for  citadels  in  the  same  format  as  in  Appendix  D  and  E  for  the  10  percent  and  40 
percent  cases. 

Table  F-l.  AREA  DEFENSE  =  100%  OF  VALUE 


Number  of  Attackers 


Value 

No  ipp 

ipp 

Difference 

Percentage 

Difference 

4,450 

11,100 

24,805 

13,750 

123.87% 

(1/4) 

8,900 

22,600 

50,800 

28,200 

124.78% 

(1/2) 

13,350 

34,250 

84,950 

50,700 

148.03% 

(3/4) 

17,800  47,050 

jm _ 
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Table  F-2.  TARGET  SET  OF  30  CITADELS 
Area  Defense  =  100  percent  of  value 


Attackers 

no  ipp 

ipp 

1000 

372 

200 

2000 

784 

400 

3000 

1214 

604.5122 

4000 

1591 

817.4713 

5000 

2042 

1017.471 

6000 

2425 

1217.471 

7000 

2836 

1393.939 

8000 

3230 

1566.130 

9000 

3622 

1737.484 

10000 

4023 

1909.048 

11000 

4412 

2080.798 

12000 

4813 

2252.357 

13000 

5202 

2424.009 

14000 

5586 

2595.466 

15000 

5992 

2868.117 

16000 

6377 

2938.561 

17000 

6756 

3110.201 

18000 

7151 

3281.761 

19000 

7534 

3453.333 

20000 

7917 

3624.869 

21000 

8299 

3796.405 

22000 

8681 

3967.983 

23000 

9064 

4139.537 

24000 

9447 

4311.188 

25000 

9830 

4482.627 

26000 

10213 

4654.273 

27000 

10595 

4825.741 

28000 

10978 

4997.365 

29000 

11361 

5168.940 

30000 

11741 

5340.476 

31000 

12124 

5512.001 

32000 

12507 

5683.561 

33000 

12890 

5855.146 

34000 

13273 

6026.652 

35000 

13645 

6198.212 

Table  F-2  (Continued) 


Attackers 

naipp 

ipp 

36000 

14028 

6369.749 

37000 

14411 

6541.273 

38000 

14781 

6712.849 

39000 

15150 

6884.418 

40000 

15536 

7055.943 

41000 

15913 

7227.484 

42000 

16280 

7399.020 

43000 

16658 

7570.545 

44000 

17041 

7742.121 

45000 

17349 

7913.646 

46000 

17600 

8058.171 

47000 

17795 

8256.610 

48000 

17800 

8428.137 

49000 

17800 

8599.755 

50000 

17800 

8771.280 

51000 

17800 

8942.657 

52000 

17800 

9114.442 

53000 

17800 

9285.490 

54000 

17800 

9455.623 

55000 

17800 

9623.533 

56000 

17800 

9790.361 

57000 

17800 

9947.421 

58000 

17800 

10109.54 

59000 

17800 

10268.95 

60000 

17800 

10426.43 

61000 

17800 

10587.74 

62000 

17800 

10740.23 

63000 

17800 

10897.48 

64000 

17800 

11050.45 

65000 

17800 

11207.02 

66000 

17800 

11356.84 

67000 

17800 

11503.34 

68000 

17800 

11649.84 

69000 

17800 

11791.78 

70000 

17800 

11926.73 

Continued 
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Table  F-2  (Continued) 


Attackers 

no  ipp 

71000 

17800 

72000 

17800 

73000 

17800 

74000 

17800 

75000 

17800 

76000 

17800 

77000 

17800 

78000 

17800 

79000 

17800 

80000 

17800 

81000 

17800 

82000 

17800 

83000 

17800 

84000 

17800 

85000 

17800 

86000 

17800 

87000 

17800 

88000 

17800 

89000 

17800 

90000 

17800 

91000 

17800 

92000 

17800 

93000 

17800 

94000 

17800 

95000 

17800 

96000 

17800 

97000 

17800 

98000 

17800 

99000 

17800 

100000 

17800 

icc 

12055.48 

12177.28 

12293.63 

12404.39 

12510.41 

12611.87 

12708.53 

12801.39 
12890.34 

12975.64 

13057.68 
13136.09 
13212.00 

13284.88 
13355.02 

13422.77 

13488.64 

13554.51 
13620.38 
13684.86 

13747.02 

13807.03 

13864.74 

13920.50 

13974.53 

14026.69 

14077.42 

14126.52 
14174.13 

14220.42 
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